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I. INTRODUCTION 

Within the standard model (SM), two-jet (dijet) events 
are produced in proton-antiproton (pp) collisions pre- 
dominantly from hard quantum chromodynamics (QCD) 
interactions of two partons. The fragmentation and 
hadronization of the outgoing partons produce hadronic 
jets. The dijet mass spectrum predicted by QCD falls 
smoothly and steeply with increasing dijet mass. Many 
extensions of the SM predict the existence of new massive 
particles that decay into two energetic partons (quarks, 
q, or gluons, <?), which can potentially be observed as a 
narrow resonance in the dijet mass spectrum. Such par- 
ticles include the excited quark axigluon [2j, flavor- 
universal coloron [|[, color-octet techni-p [3], Randall- 
Sundrum (RS) graviton [f|, W, Z' @, and diquark in 
the string-inspired Eq model [7j. 

Here we briefly discuss the theoretical models for these 
new particles. In the SM, the quarks are considered as 
fundamental particles. However, the presence of their 
generational structure and mass hierarchy motivate mod- 
els of quark compositeness in which the quarks consist 
of more fundamental particles. If a quark is a compos- 
ite particle, an excited state of a quark q* is expected, 
which decays to qg In chiral color models, the SU(3) 
gauge group of QCD results from the spontaneous break- 
ing of the chiral color gauge group of SU(3) x SU(3). 
Any model of chiral color predicts the presence of the 
axigluon, a massive axial vector gluon, that decays to 
qq [2| . The flavor- universal coloron model also embeds 
the SU(3) of QCD in a larger gauge group and predicts 
the presence of a color-octet coloron which decays to 
qq Q. Technicolor models seek to explain electroweak 
symmetry breaking via the dynamics of new interactions 
among techniquarks. The models of extended techni- 
color and topcolor-assisted technicolor predict the pres- 
ence of a color-octet techni-p (pra) which decays to qq 
or gg Q . The RS model of a warped extra dimension of- 
fers a solution for the hierarchy between the electroweak 
scale and Planck scale Mpi by introducing an extra spa- 
cial dimension This model predicts a Kaluza-Klein 
tower of graviton states (RS gravitons) which decay to 
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qq or gg. The grand unified theories (GUT) based on 
larger gauge groups, e.g., E$ and 50(10), or left-right- 
symmetric models jgj often introduce additional gauge 
bosons, such as W and Z' , which decay to qq' and qq, 
respectively [6]. The E 6 GUT model also predicts the 
presence of a diquark which decays to qq or qq 0- 

In the past, the UA2 [9], CDF 0, 03 and DO [13 
experiments searched for resonances in the dijet mass 
spectrum and set limits on their production. In this ar- 
ticle, we present a first measurement of the dijet mass 
spectrum and a search for massive particles which decay 
into dijets in pp collisions at the center-of-mass energy 
\/i =1.96 TeV. This analysis uses data corresponding to 
an integrated luminosity of 1.13 fb _1 collected between 
February 2002 and February 2006 with the CDF II de- 
tector at the Fermilab Tevatron. 

The measurement of the dijet mass spectrum is also an 
important test of perturbative QCD (pQCD) predictions. 
It provides complementary information to the inclusive 
jet cross section measurements [H, [U, [IH, [ljl [13], and 
comparisons of the measurement with pQCD predictions 
provide constraints on the parton distribution functions 
(PDFs) of the proton, in particular at high momentum 
fraction x (x > 0.3) where the gluon distribution is not 
well constrained 11811. 



II. THE CDF DETECTOR 

The CDF II detector is described in detail else- 
where 19]. Here, the components that are relevant 
to this search are briefly described. Surrounding the 
beam pipe, there is a tracking system consisting of a 
silicon microstrip detector, a cylindrical drift chamber, 
and a solenoid magnet that provides a 1.4 T magnetic 
field. The central and plug calorimeters, which cover 
the pseudorapidity regions of jr?] < 1.1 and 1.1 < 
|ry| < 3.6 (20l |. respectively, surround the tracking sys- 
tem with a projective tower geometry and measure the 
energy of interacting particles. The calorimeters are 
segmented into electromagnetic and hadronic sections 
that consist of lead-scintillator and iron-scintillator, re- 
spectively. In the central region, the calorimeter con- 
sists of 48 modules, segmented into towers of granular- 
ity At] x A<f> w 0.1 X 0.26. The energy resolution of 
the central electromagnetic calorimeter for electrons is 
a(E T )/E T = 13.5%/^E T (GeV)ffil.5%, while the energy 
resolution of the central hadron calorimeter for charged 
pions that do not inter act in the electromagnetic section 
is <j{E t )/E t = W%U/E T {GeV) © 3%, where E T is the 
transverse energy [201 ] . The wall hadron calorimeter cov- 
ers the gap in the projective tower geometry between the 
central and plug hadron calorimeters, corresponding to 
0.7 < \rj\ < 1.3, with segmentation similar to that of the 
central calorimeter. The energy resolution of the wall 
hadron calorimeter is <t(E t )/E t = 75%/ v / E T (GeV) © 
4% for charged pions that do not interact in the electro- 
magnetic section. A system of Cherenkov counters, lo- 
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cated around the beam pipe and inside the plug calorime- 
ters, is used to measure the number of inelastic pp colli- 
sions per bunch crossing and thereby the luminosity. 



III. DATASETS AND EVENT SELECTION 

The dijet data used in this search were collected us- 
ing a three-level online event selection (trigger) system, 
and are identical to the data used in Ref. [17|. The 
four trigger selections used in this analysis are referred 
to as "jet20", "jet50", "jet 70", and "jet 100" according 
to the Et threshold (in GeV) of calorimeter clusters re- 
constructed using a cone algorithm 21| with cone ra- 
dius R cone = 0.7. The jet20, jet50, and jet70 trigger 
rates are randomly reduced (prescaled) to avoid saturat- 
ing the bandwidth of the data acquisition system. Af- 
ter prescaling, these trigger datasets correspond to 1.44, 
32.5, and 143 pb _1 of integrated luminosity, respectively. 
The jetlOO trigger is not prescaled. 

Jets are reconstructed from the energy depositions in 
the calorimeter towers with the transverse momentum 
Pt above 0.1 GeV/c. Jets are formed from the four- 
vectors of calorimeter towers [22j using the cone-based 
midpoint jet clustering algorithm [lj [23, [24[ with cone 
radius R CO ne = 0.7. The kinematics of a jet is defined by 
the four- vector recombination scheme [23j. 

Cosmic ray and beam lo ss back ground events are re- 
moved by requiring pr / &t < min(3 + 0.0125 x 
p^f* ,6), where Ifir and ^2 Et are the missing Et and 
sum Et HH, respectively, and p^ 1 (GeV/c) is the pt 
of the leading jet in the event before the corrections de- 
scribed in Sec. IIVI are applied. This requirement makes 
the fraction of background events negligible, and has an 
efficiency of > 95% for dijet events. In order to ensure 
good coverage of each event by the detectors, the primary 
event vertex is required to be within 60 cm of the cen- 
ter of the detector along the z-axis [2(|. The efficiency 
for this requirement is determined to be 96% from the 
distribution of primary event vertices along z measured 
using a sample of minimum bias events. 



IV. JET CORRECTIONS AND 
MEASUREMENT OF THE DIJET MASS 
SPECTRUM 

The jet energies measured by the calorimeters are af- 
fected by instrumental effects such as calorimeter non- 
uniformity, non-linearity, and energy smearing. We cor- 
rect for these biases in several steps (2(|. First, an 77- 
dependent relative correction is applied to equalize the 
response of the calorimeter. The equalized jet pt is 
then corrected for the effects of pileup, i.e., additional 
pp interactions in the same bunch crossing. The pileup 
correction subtracts 0.97 ± 0.29 GeV/ c for each addi- 
tional primary vertex from the measured jet pt- Then, 
a pT-dependent correction is applied to account for, on 



average, the under-measured hadron energy due to the 
non-linearity of the calorimeter response. The correction 
factors are 1.19 and 1.06 at jet pt — 90 and 600 GeV/c, 
respectively. After these corrections we reconstruct the 
dijet mass rrijj from the four-vectors of the two highest 
corrected-pT jets. We form the dijet mass spectrum as: 



da 



dm 



j j 



1 



(i) 



where m is the observed number of events, ttrig,i is the 
trigger efficiency, Amjj_i is the bin width, and Ci is the 
integrated luminosity of the trigger dataset used for the 
i-th dijet mass bin. The bin width is set to 10% of the 
dijet mass which approximately corresponds to the dijet 
mass resolution. We count only events in which both of 
the leading two jets have a rapidity magnitude \y\ (20j 
less than 1. We use the jet20, jet50, jet70, and jetlOO 
datasets for the dijet mass regions of 180-241, 241-321, 
321-427, and above 427 GeV/c 2 , respectively, where the 
trigger efficiencies are higher than 99.8%. This ensures 
a negligible uncertainty from trigger efficiency measure- 
ments. 

In order to compare the measured dijet mass spectrum 
with QCD predictions at the hadron level [13], the mea- 
sured spectrum must be further corrected for the bin-by- 
bin migration effect due to the finite resolution of the 
rrijj measurement and the efficiencies of the offline event 
selection requirements. We obtain this correction using 
QCD dijet events generated by the PYTHIA 6.2 [23] Monte 
Carlo simulation program, that have passed through the 
CDF detector simulation (28|. The pythia events are 
generated with Tune A [2§] , which refers to the set of pa- 
rameters describing multiple-parton interactions and ini- 
tial state radiation that have been tuned to reproduce the 
energy observed in the region transverse to the leading 
jet [30(. It has also been shown to provide a reasonable 
description of the measured energy distribution inside a 
jet [311 ]. The correction is determined on a bin- by-bin ba- 
sis by taking the ratio of a hadron-level cross section to 
a calorimeter-level cross section. The hadron-level cross 
section is defined using hadron-level jets clustered from 
the final state stable particles [13] in PYTHIA with the 
same jet clustering algorithm as the one used to cluster 
calorimeter towers. The leading two jets are required to 
have \y\ < 1. The calorimeter-level cross section is ob- 
tained by analyzing the pythia events using the same 
analysis chain as for the data. Since the correction de- 
pends on the dijet mass spectrum, the pythia events are 
reweighted to match the dijet mass spectrum measured 
in data before the correction factor is calculated. The 
size of the correction ranges from ~ 1.2 at low rrijj (200 
GeV/c 2 ) to - 1.5 at high m 3] (1250 GeV/c 2 ). 

The systematic uncertainties arise mainly from four 
sources: the jet energy scale, the jet energy resolution, 
the unfolding correction, and the integrated luminosity. 
The dominant source is from the absolute jet energy 
scale. The size of the uncertainty in the cross section 

% at high rrijj. The 



varies from 10% at low rrijj to 
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FIG. 1: (a) The measured dijet mass spectrum for both jets 
to have \y\ < 1 compared to the NLO pQCD prediction ob- 
tained using the CTEQ6.1 PDFs. (b) The ratio of the data 
to the NLO pQCD prediction. The experimental systematic 
uncertainties, theoretical uncertainties from PDF, the ratio of 
MRST2004/CTEQ6.1, and the dependence on the choice of 
renormalization and factorization scales are also shown. An 
additional 6% uncertainty in the determination of the lumi- 
nosity is not shown. 



uncertainty in the relative jet energy scale introduces a 
3% uncertainty on the measured cross section at low rrijj 
and lg % at high rrijj. The uncertainty in the modeling 
of jet energy smearing is estimated from the difference 
between the data and PYTHIA samples using the bisector 
method [H[ , and it introduces an uncertainty in the cross 
section of 1% at low rrijj and +5% at high rrijj. The dif- 
ference between the unfolding correction from a pythia 
sample and that from a sample generated by the herwig 
6.5 [3J] Monte Carlo simulation program is also taken as 
a systematic uncertainty to account for the uncertainty 
in the modeling of jet fragmentation. This uncertainty is 
2% at low rrijj and 8% at high rrijj. The uncertainty in 
the determination of the integrated luminosity is 6%, in- 
dependent of rrijj. The total uncertainty is at low 



j j 



and _ 49 % at high rrijj. The measured dijet mass 
spectrum after all the corrections discussed above are 
applied is shown together with the statistical and sys- 
tematic uncertainties in Fig.[TJa). 



V. COMPARISONS WITH QCD PREDICTIONS 

The measured dijet mass spectrum is compared in 
Fig. Q] to the next-to-leading-order perturbative QCD 
(NLO pQCD) predictions from fastNLO H. The Pre- 
dictions were obtained using the CTEQ6.1 [lj§] PDFs 
with the renormalization and factorization scales both 
set to /io, the average pr of the leading two jets. Jets 
are reconstructed by the midpoint algorithm with cone 
radius R CO ne = 0.7. The maximum separation between 
two partons merged into a jet is set to R CO ne x Rsep [36T ]. 
and R S ep is set to 1.3 



211 . I24j . Setting the renormal- 
ization and factorization scales to 2/xrj instead of fio re- 
duces the cross section prediction by 5-10%, and setting 
R sep — 2 increases the cross section by < 10%. The PDF 
uncertainties estimated from 40 CTEQ6.1 error PDFs 
and the ratio of the predictions using MRST2004 37J and 
CTEQ6.1 are shown in Fig.[TYb). The PDF uncertainty 
is the dominant theoretical uncertainty for most of the 



rrijj range. The NLO pQCD predictions for jets clustered 
from partons need to be corrected for non-perturbative 
underlying event and hadronization effects. The multi- 
plicative parton-to-hadron- level correction (C p ^h) is de- 
termined on a bin-by-bin basis from a ratio of two dijet 
mass spectra. The numerator is the nominal hadron-level 
dijet mass spectrum from the pythia Tune A samples, 
and the denominator is the dijet mass spectrum obtained 
from jets formed from partons before hadronization in 
a sample simulated with an underlying event turned off. 
We assign the difference between the corrections obtained 
using herwig and pythia Tune A as the uncertainty on 
the C p ^h correction. The C p -,h correction is 1.16 ±0.08 
at low rrijj and 1.02±0.02 at high rrijj. Figurc[T]shows the 
ratio of the measured spectrum to the NLO pQCD pre- 
dictions corrected for the non-perturbative effects. The 
data and theoretical predictions are found to be in good 
agreement. To quantify the agreement, we performed a 
X 2 test which is the same as the one used in the inclu- 
sive jet cross section measurements [HI, The test 
treats the systematic uncertainties from different sources 
and uncertainties on C p —>h as independent but fully cor- 
related over all rrijj bins and yields x 2 / n -d.f. = 21/21. 



VI. SEARCH FOR DIJET MASS RESONANCES 

We search for narrow mass resonances in the measured 
dijet mass spectrum by fitting the measured spectrum to 
a smooth functional form and by looking for data points 
that show significant excess from the fit. We fit the mea- 
sured dijet mass spectrum before the bin-by-bin unfold- 
ing correction is applied. We use the following functional 
form: 



da 
dm 



Po (l-x) Pl /x P2+P3Mx \ x = m jj /Vs, (2) 



J3 



where po? Pi, P2, and are free parameters. This form 
fits well the dijet mass spectra from pythia, herwig, 
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FIG. 2: (a) The measured dijet mass spectrum (points) fitted 
to Eq. ([2]) (dashed curve). The bin-by-bin unfolding correc- 
tions is not applied. Also shown are the predictions from the 
excited quark, q* , simulations for masses of 300, 500, 700, 
900, and 1100 GeV/c 2 , respectively (solid curves), (b) The 
fractional difference between the measured dijet mass distri- 
bution and the fit (points) compared to the predictions for 
q* signals divided by the fit to the measured dijet mass spec- 
trum (curves). The inset shows the expanded view in which 
the vertical scale is restricted to ±0.04. 



and NLO pQCD predictions. The result of the fit to 
the measured dijet mass spectrum is shown in Fig. [21 
Equation (T5|) fits the measured dijet mass spectrum well 
with % 2 /n.d.f. = 16/17. We find no evidence for the 
existence of a resonant structure, and in the next section 
we use the data to set limits on new particle production. 



VII. LIMITS ON NEW PARTICLE 
PRODUCTION 

Several theoretical models which predict the existence 
of new particles that decay into dijets are considered in 
this search. For the excited quark q* which decays to qg, 
we set its couplings to the SM SU(2\ U(l), and SU(3) 
gauge groups to be / = /' = f s = 1 [l[ , respectively, and 
the compositeness scale to the mass of q*. For the RS 
graviton G* that decays into qq or gg, we use the model 
parameter k/Mpi =0.1 which determines the couplings 
of the graviton to the SM particles. The production cross 
section increases with increasing k/Mpi; however, values 
of k/Mpi > 0.1 are disfavored theoretically [11]. For W 
and Z' , which decay to qq' and qq respectively, we use 
the SM couplings. The leading-order production cross 
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FIG. 3: Dijet mass distributions for simulated signals of the 
q*, RS graviton, W', and Z' with the mass of 800 GeV/c 2 . 



sections of the RS graviton, W', and Z' are multiplied by 
a factor of 1.3 to account for hig her-order effects in the 
strong coupling constant a s [391] . All these models are 
simulated with pythia Tune A. Signal events of these 
models from pythia are then passed through the CDF 
detector simulation. 

The dijet mass distributions from q* simulations with 
masses 300, 500, 700, 900, and 1100 GeV/c 2 are shown 
in Fig. [21 The dijet mass distributions for the q* , RS 
graviton, W' , and Z' simulations with the mass of 800 
GeV/c 2 are shown together in Fig. [3l The shapes of the 
distributions are mainly determined by the jet energy 
resolution and QCD radiation which leads to tails on the 
low mass side. Since the natural width of these parti- 
cles is substantially smaller than the width from the jet 
energy resolution, all the dijet mass distributions appear 
similar. However, the dijet mass resonance distributions 
are somewhat broader for q* and RS gravitons than for 
W' and Z' because q* and RS gravitons can decay into 
the mode containing gluons, unlike W and Z' . Gluons 
radiate more than quarks and tend to make the resulting 
dijet mass distributions broader. As a result, the cross 
section limits obtained based on the q* and RS gravi- 
ton resonance shapes are about 20% larger than those 
obtained with the W and Z' resonance shapes. 

We also consider production of the axigluon A that 
decays into qq, E 6 diquark D (D c ) that decays into 
qq (qq), and color-octet techni-p (prs) that decays into 
qq or gg. Their lowest-order theoretical predictions for 
fjsig = (x ■ B ■ A are shown in Fig. 0] along with the pre- 
dictions for the other models described above as a func- 
tion of new particle mass, where a is the new particle 
production cross section, B is the branching fraction to 
dijets, and A is the kinematical acceptance for each of 
the leading two jets to have \y\ < 1. In addition, the 
flavor-universal coloron C which decays to qq is consid- 
ered. The cross section for the coloron is always larger 
than or equal to that for the axigluon, so the limits on 
the axigluon apply to the coloron as well. For pxs pro- 
duction, predictions are for the mass-degenerate pxa with 
the standard topcolor-assisted-technicolor couplings and 
with the set of parameters in [iol ]. 

We set upper limits on a slg = cr-B-A as follows. We use 
the likelihood function L — Y\ i ex P( — A*i)/ n i'j where 
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TABLE I: Observed 95% C.L. upper limits on the new particle 
production cross sections times branching fraction to dijets 
times the acceptance for both jets to have \y\ < 1 obtained 
with the signal shapes from W' , Z' , RS graviton (G*), and 
q* production. 



Mass 
(GeV/c 2 ) 



260 
280 
300 
320 
340 
360 
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400 
425 
460 
500 
540 
580 
620 
660 
700 
750 
800 
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900 
950 
1000 
1050 
1100 
1150 
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1250 
1300 
1350 
1400 



W 
1.1 x 10 2 

8.1 x 10 1 

4.5 x 10 1 
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G* 



1.1 x 
8.3 x 
5.1 x 
3.1 x 
1.9 x 
1.1 x 
8.6 x 
7.3 x 
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Hi = a sl9 Citing 9 /n^ l 9 +nf D is the predicted number of 
events, Ci is the event selection efficiency in the i th dijet 
mass bin, and 9 /n* f is the predicted fraction of signal 
events in bin i. We model the QCD dijet mass spectrum 
with Eq. ([2]) and use Eq. ([1]) to extract nf CD from the 
differential cross section da/dmjj. For each value of a slg 
we maximize the likelihood with respect to the four pa- 
rameters in Eq. @. We integrate this profiled likelihood 
over Bayesian priors for the parameters describing the 
systematic uncertainties [4l|, and we use a flat prior on 
a sl9 to extract Bayesian upper limits on that parame- 
ter. Although this procedure uses Bayesian techniques, 
we verified that the resulting upper limits have good fre- 
quentist coverage. 

The obtained 95% confidence level (C.L.) limits us- 
ing the W, Z', RS graviton, and q* signal resonance 
shapes are shown in Fig. [4] and Table [J as a function 
of the new particle mass. Also shown in Fig. 0] are the 
theoretical predictions for the various models. For the 
W, Z', q* , and RS graviton, the mass exclusion is de- 
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FIG. 4: Observed 95% C.L. upper limits on new particle pro- 
duction cross sections times the branching fraction to dijets 
obtained with the signal shapes from (a) W' , (b) Z' , (c) RS 
graviton, and (d) q* production. Also shown are the cross sec- 
tion predictions for the production of W 1 , Z\ RS graviton, 
Pts, 9*, axigluon, flavor-universal coloron, and diquark 
for the set of parameters described in the text. The limits 
and theoretical predictions are for events in which both of 
the leading two jets have \y\ < 1. 



termined by comparing the limits obtained with their 
respective signal shapes and a slg predictions. For the 
axigluon, flavor-universal coloron, and Eq diquark, we 
compare their a sl9 predictions to the limits obtained with 
the q* signal shapes; these particles do not decay into the 
mode containing a gluon, so their signal shape would be 
narrower than that of q* , and thus the mass exclusions 
obtained with the q* signal shapes are conservative. For 
PT8i w e compare its a s%9 predictions with the limits ob- 
tained using the RS graviton signal shapes. The prs and 
RS graviton decay channels are similar; the branching 
fraction to the gg state is higher for the RS gravitons 
than for pxs , so comparing the limits obtained using the 
RS graviton signal shapes to the pts's <J slg predictions 
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yields a conservative mass exclusion. 

The mass exclusion regions obtained in this search are 
260 < m < 870 GeV/c 2 for g*, 260 < m < 1250 GeV/c 2 
for the axigluon and flavor-universal coloron, 290 < m < 
630 GeV/c 2 for the E 6 diquark, 260 < m < 1100 GeV/c 2 
for pTs, 280 < m < 840 GeV/c 2 for W, and 320 < m < 
740 GeV/c 2 for Z'. For the RS graviton, this search did 
not exclude any mass region. As the mass exclusions from 
the previous dijet mass resonance searches [13, [HI G3 
are m < 775 GeV/c 2 for q*, m < 980 GeV/c 2 for the 
axigluon and flavor-universal coloron, 290 < m < 420 
GeV/c 2 for the E e diquark, and 260 < m < 480 GeV/c 2 
for pts, this search provides the most stringent lower 
mass limits for these particles. This search also extends 
the existing mass exclusions of 300 < m < 800 GeV/c 2 
for W and 400 < m < 640 GeV/c 2 for Z' obtained in the 
previous dijet mass resonance search (l2| . However, the 
more restrictive mass exclusions of m < 1000 GeV/c 2 for 
W and m < 923 GeV/c 2 for Z' come from the W -> ezv 
search [42j] and Z' — > e + e _ search [43| . 

VIII. CONCLUSIONS 

In conclusion, we performed a search for new particles 
which decay into dijets by measuring the dijet mass spec- 
trum using pp collision data from 1.13 fb _1 of integrated 
luminosity collected with the CDF II detector. The mea- 
sured dijet mass spectrum is found to be consistent with 
NLO pQCD predictions based on recent PDFs and does 
not show evidence of a mass resonance from new particle 
production. We set 95% confidence level upper limits on 



new particle production cross sections times the branch- 
ing fraction to dijets times the acceptance for both jets 
to have \y\ < 1. We also determine the mass exclusions 
for the excited quark, axigluon, flavor-universal coloron, 
Eq diquark, color-octet techni-p, W', and Z' for a spe- 
cific representative set of model parameters. This search 
sets the most stringent lower mass limits on the excited 
quark, the axigluon, the flavor-universal coloron, the Eq 
diquark, and the color-octet techni-p. 
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